function y = func_p_from_q(q, CONSTS)
    eps_a = CONSTS.eps_a;
    %y = -1i*q; 
    y = sqrt_from_negative_num(eps_a-q.^2); %(1-q.^2).^0.5;
    
end
